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Abstract 

We present a definition for second-order freeness in the quaternionic 
case. We demonstrate that this dehnition on a second-order proba¬ 
bility space is asymptotically satished by independent symplectically 
invariant quaternionic matrices. 

This definition is different from the natural definition for com¬ 
plex and real second-order probability spaces, those motivated by the 
asymptotic behaviour of unitarily invariant and orthogonally invariant 
random matrices respectively. 

Most notably, because the quaternionic trace does not have the 
cyclic property of a trace over a commutative field, the asymmetries 
which appear in the multi-matrix context result in an asymmetric con¬ 
tribution from the terms which appear symmetrically in the complex 
and real cases. 


1 Introduction 


In 24 , D.-V. Voiculescu established a connection between free probability, a 


noncommutative analogue of classical probability, and the behaviour of large 
random matrices. Many ensembles of large random matrices exhibit this 
asymptotic freeness, a property characterizing the behaviour of the expected 
value of the trace of a product of independent matrices in terms of the 
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behaviour of the matrices individually. The same behaviour is observed in 
complex, real, and quaternionic random matrices. 

Second-order freeness is developed in [7,14 15 in order to similarly char¬ 
acterize the covariances of traces of multi-matrix expressions in high dimen¬ 
sions. The behaviour described here is exhibited by independent ensembles 
of matrices whose probability distribution is unitarily invariant (unchanged 
if the matrices are conjugated by an arbitrary unitary matrix), a natural 
symmetry for complex random matrices. The behaviour of orthogonally in¬ 
variant random matrices (a natural symmetry for real matrices) is discussed 
in [21| . Independent orthogonally invariant random matrices do not gen¬ 
erally satisfy the definition given in 15 , so another definition is provided 


for real second-order freeness. This paper gives a definition of quaternionic 
second-order freeness, motivated by the covariance of traces of products of 
large symplectically invariant matrices. 

We consider the covariance of the quaternionic trace (as opposed to the 
real part of the trace). While the behaviour of quaternionic eigenvalues are 
more complicated than eigenvalues over a commutative field (for example, 
the qualitatively different behaviour of left and right eigenvalues), eigen¬ 
values may nonetheless be quaternion-valued, and the sum of quaternionic 
eigenvalues is related to the sum of the quaternionic diagonal entries rather 
than to the real part of this sum (see, e.g., [^, Chapter 9). Thus the 
quaternion-valued (partial) trace is the relevant quantity. 

If E(Tr(Xfc)) = E(Tr(17;)) = 0, the asymptotic covariance will vanish 
if the number of independent terms in the two traces is different, and if the 
number is the same: 


lim cov (Tr {Xi ■ ■ ■ Xp ), Tr (W • • • Yp)) 
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where subscripts are taken modulo p (see Definition 2.31 for the full def¬ 
inition). This is different from the analogous expression for unitarily and 
orthogonally invariant matrices in two notable ways. 

Most notably, because the quaternionic trace is not cyclic on products 
(because the quaternions are not commutative), the cyclic symmetries of 
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traces break down in the multi-matrix context which is integral to the free¬ 
ness and second-order freeness setting. The real and complex expressions 
are sums over terms known as “spoke diagrams” (see Figure 2.6), which 
are invariant under cycling the terms within each of the two traces. In 
the quaternionic case, two distinguished terms appear with a different co¬ 
efficient. (In the real case, all terms appear with the coefficient 1. In the 
complex definition, the first two terms appear with coefficient 1, and the 
last two terms do not appear.) 

Secondly, the contribution from each “spoke” is the expected value of 
the real part of a trace, rather than simply the expected value of the trace. 
The real-part function does not appear in either the real or the complex 
case, where the covariance may be complex-valued when a unitarily invari¬ 
ant or complex-valued orthogonally invariant matrix appears (for example 
a Wishart matrix (see Definition 5.4) G*DG or G^DG where the determin¬ 
istic matrix D is complex-valued). However, quaternionic conjugation is a 
diffeomorphism (unlike complex conjugation), so the contribution of a term 
and its conjugate in a symplectically symmetric distribution will be equal, 
and their quaternion-imaginary parts will cancel. 

In this paper, we use a topological expansion for quaternionic multi¬ 
matrix expressions from (see also [^[^). Asymptotic terms are char¬ 
acterized by the annular noncrossing conditions in [13] (see also M 21 for 


more on the unoriented case). The necessary background is outlined in Sec¬ 
tion The definitions of quaternionic second-order freeness and asymptotic 
quaternionic second-order freeness are also given in this section. The asymp¬ 
totic behaviour of cumulants of traces is given in Section The proof that 
independent symplectically invariant random matrices satisfy the definition 
of quaternionic second-order freeness is given in Section A number of 
important random matrix ensembles (Ginibre, GSE, Wishart, and Haar- 
distributed symplectic) are discussed in Section including demonstrations 
that they satisfy the necessary convergence conditions and some information 
on the values of their contributions to expression Q. 


2 Preliminaries 

2.1 Partitions, permutations, and maps 

Definition 2.1. For m,n G N, we let [n] := {!,...,n} and [m,n] = 
{m,... ,n}. 

We will sometimes append a formal point at infinity to a set: loo = 
I U {oo}. 
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For a set / C Z, we let —I := {—k : /c G /} and ±1 := / U {—I). 

Definition 2.2. A partition (or a set partition) tt = {Vi, ..., Vm} of a set 
/ is a set of nonempty subsets of I (called blocks of the partition) such that 
Vi n = 0 for all i ^ j, and IJi^i Vi = I- We denote the set of partitions of 
I hy V (!) and the set of partitions of [n] by V (n). 

For TT = {Vi,...,Vm} ev (/), we let ivr := {iVd,. . .,±Vm} G V (±1). 
If J C /, we let vr| j = {Vi n J, ..., Vm Ci J} G V (J) (ignoring empty blocks). 
We denote the number of blocks of partition tt G P (/) by ^ (tt). 

For a function /;/—)■ J, we dehne the kernel of f, denoted ker (/), to 
be the partition of I whose blocks are the nonempty preimages f~^ (j) for 
each j G J. 

The partitions V (I) form a poset, where, for tt,p G V (I), we say that 
TT ^ p if every block of vr is a subset of a block of p. 

The smallest element of the poset V (I) (in which each block contains 
exactly one element) is denoted 0/, and the largest element (in which all the 
elements are contained in one block) is denoted 1/. We denote 0^ := 0[„] 
and In := l[n]) or simply 0 and 1 if / is understood. 

The join of two partitions TT,p G V (/), denoted tt V p, is the smallest 
element which is bigger than both vr and p. The meet of vr and p, denoted 
vrAp, is the largest partition which is smaller than both tt and p. (The poset 
of partitions forms a lattice, so the meet and join are uniquely defined.) 

We say that a partition tt G V (I) connects blocks V,W G p of another 
partition p gV {!) \iV and W are subsets of the same block of tt V p. 

A partition whose blocks all contain exactly 2 elements is called a pairing. 
We denote the set of all pairings of I by 1^2 {!), and the set of all pairings 
of [n] by V 2 (n). 

The following lemma will be useful: 

Lemma 2.3. If tt, p, a G V (/), then {tt y p) — {tt y p y a) < # (vr) — 
#(7r Vcr). 

Proof. If k blocks of tt V p are joined by a block of a, then at least as many 
blocks of TT are joined by that block. Thus joining with a block of a reduces 
the number of blocks of tt by at least much as it reduces the number of 
blocks of TT V p. Repeating for all the blocks of a, the result follows. □ 

Definition 2.4. We define an integer partition A of integer n > 0 as a list 
of integers (parts) Ai > • • • > > 0 such that Ai + • • • + A^ = n. 
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If it is unclear from the context, we will specify whether we are referring 
to a set partition or an integer partition. If it is not stated, a partition is a 
set partition. 


Definition 2.5. We denote the set of permutations on set / by S' (I) and 
the permutations on [n] by S'„. 

We will write permutations in cycle notation (sometimes omitting cycles 
with only one element). We follow the convention that permutations are 
evaluated right-to-left. 

We note that conjugation by another permutation substitutes each ele¬ 
ment in the cycle notation with its image under the conjugating permutation, 
and thus that conjugation by another permutation preserves the cycle struc¬ 
ture of the permutation (and thus cycling the permutations in a product also 
preserves the cycle structure). 

The orbits, or cycles, of a permutation form a partition of the domain 
set. We will also write the number of cycles of a permutation vr (that is, the 
number of blocks in this partition) as ^ (vr). The domain will be specified 
if it is not clear from the context. We will also use other partition notation, 
such as the join V, implicitly applied to the orbits of permutations. (This 
may be ambiguous in the definition of PMnc (vr). Definition 2.10 While in 
PM (tt), tt is interpreted as a partition, in PM^c (vr) it is a permutation.) 

We will also consider tt € V 2 (I) to be permutations, where an element 
is taken to the element it is paired with. 

If TT G S' (I) and J C /, we dehne the permutation induced by vr on 
J, which we denote by 7r| j, by letting 7r| j (k) = vr™, where m > 0 is the 
smallest integer such that vr™ G J. In cycle notation, this amounts to 
deleting elements which are not in J. If none of the cycles of vr contain 
elements of both J and I \ J, then this is simply the restriction of the 
function to J. 

We will frequently subscript functions by a permutation. Let / be a 
function which is invariant under cycling of terms in a product (such as 
the trace of matrices over a commutative field). Let I C ± [n], and let 
TT = (ci,... ,Cni) • • • (cn^_i, • • •, Cn^) £ -S' (/). Then 


fn (xi, . . . ,Xn) := / (5 


'Cl 


‘Xc 


..)■■■/( 




Multiplying (left or right) a permutation tt by a transposition (a permu¬ 
tation of the form (a, h)) increases (vr) by 1 if a and b connect two cycles 
of TT, and decreases # (tt) by 1 if they do not. 

We use permutations to encode graphs in surfaces (see Example 2.8). 
Faces, edges, or vertices are encoded as a permutation by enumerating the 
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ends of the edges (or hyperedges) which appear around each face (or hyper 
edge or vertex) in cyclic order in a cycle of the permutation. See [9 12 
more detail on the construction, (l9 22 for more on the construction in the 


for 


unoriented case, and [1,17,20 for more on the quaternionic case. 

For any surface, we may construct an orientable two-sheeted covering 
space by letting points be close in the covering space if they are close in the 
base space by an orientation preserving path. For more on this construction, 
see [^, Section 3.3. We use k and —k to represent the two preimages of k. 
Objects in the covering space should be consistent with the covering map, 
as described in the following definition: 


Definition 2.6. We will denote the function k i—?• —k by 6. 

We will denote by PM (/) the set of all permutations vr on ±1 such that 
dirS = Tr~^ and such that no cycle of vr contains both k and —k. We denote 
PM([n]) by PM(n). 

For p G 7^(7), we denote by PM (p) the set of vr G PM (!) such that 
every cycle of vr is contained in a block of ±p. 

If vr G PM(/), the cycles are paired: for each cycle there is another 
where the elements with the signs reversed appear in reverse cyclic order. 
We choose the cycle from each pair in which the element of smallest absolute 
value (or inhnity) appears as a positive integer. We define FD (vr) as the 
permutation which is the product of these cycles, taking the domain as the 
set of elements appearing in those cycles. We will also use this notation to 
refer to the domain as a set. 


It is possible to calculate the permutation representing the vertices from 
the permutations representing the faces and the hyperedges: 

Definition 2.7. If G S (I) and a G PM (I), we let (p_ := 6ip^5. We 
define 

K ((p+,a) := 

In [^ it is shown that K ((p+, a) G PM (/). 

Letting cp := we define the Euler characteristic: 

X ((P+, a) := # ((p) /2 + # (a) /2 + # (iF (p, a)) /2 - |/| . 

We note that while K ((p+, a) depends on which cycles of (p come from (p+ 
and which from (p_, the number of cycles, and hence the Euler characteristic, 
does not, so we may write y ((p, a). 
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Figure 1: A hypermap on three faces. The middle face is shown from behind. 


Example 2.8. In Figure an example of a hypermap is shown. The middle 
face is shown from behind, so the integers are negative and appear in reversed 
(clockwise) order. The faces are represented by the permutation 

(^+ = (1, 2, 3,4, 5, 6) (7, 8,9,10,11,12,13) (14,15,16). 

Permutation gives us the fronts of the faces, while permutation 

= (-6, -5, -4, -3, -2, -1) (-13, -12, -11, -10, -9, -8, -7) (-16, -15, -14) 

gives the backs of the faces. 

The hyperedges are represented by permutation 

a = (1) (-1) (2, 6, 5) (-5, -6, -2) (3, -11, -13) (13,11, -3) 

(4, -7, -8, -9, -10) (10,9,8, 7, -4) (12) (-12) (14,15) (-15, -14) (16) (-16) . 

(We use the convention that all elements (faces, hyperedges, vertices) are 
oriented counter-clockwise internally, which is slightly different from the 
conventions used in some of the cited references.) 

The integer representing the corner of a face appearing in a given vertex 
is the integer labelling the edge of which it is the counter-clockwise vertex 
(so the clockwise vertex if the face is viewed from behind). We calculate 
that the vertices will be represented by permutation 

K (7>-H,a) = 75+^a"V- 

= (1,6) (-6, -1) (2, 4, -10) (10, -4, -2) (3, -13) (13, -3) (5) (-5) (7) (-7) 

(8) (-8) (9) (-9) (11,12) (-12, -11) (14) (-14) (15,16) (-16, -15). 
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We compute that the Euler characteristic is 


X a) = # {ip) /2 + # (a) l2 + i^{K {p, a))/2-\I\ 

= 3 + 7 + 10-16 = 4. 

This is consistent with the intuition that the diagram is noncrossing, and 
thus each of the two connected components is a sphere. 

The Euler characteristic is at most 2 times the number of connected 
components: 

Lemma 2.9. Let p, a € FM (I). Then 

x{p,oi) < 2#{±pVa) . 


See 21 for a proof. 


Definition 2.10. Let p+ G S{I), with p G PM (I) defined as above. 
We denote the set of a G PM (!) such that x (++)«) = 2# (±++ V a) by 
PMnc (++) or by PMnc (+)■ 


In some contexts, we may characterize a G PMnc (+) in terms of crossing 
conditions. In particular, we may do so when ^ (++) is equal to 1 or 2. 

Definition 2.11. Let p,a € S (I) where ^ [p) = 1. 

We say that a is disc nonstandard on p if there are distinct a,b,c G I 
such that p\{a,b,c} = «l{aAc} = 

We say that a is disc crossing on p if there are distinct a,b,c,d G I such 
that p\{a,b,c,d} = {a,b,c,d) and a|{n^b,c,d} = («> c) {b,d). 

We say that a is disc noncrossing on p if it is neither disc nonstandard 
or disc crossing. We denote the set of a G 5 (I) which are disc noncrossing 
on p by S'disc-nc (+)■ 

Definition 2.12. Let p,a G S {!), where #{p) = 2. We say that a is 
annular nonstandard on p if it satisfies at least one of the following two 
conditions: 


1. There exist distinct elements a,b,c G I snch that p\^^ ^ ^ c} ~ 

(a,5,c). 

2. There exist distinct elements a,b,c,d G I such that p\^abcd} ~ (“’ b) (c, d) 

and = {a,c,b,d). 
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For X, y G / in distinct cycles of we define 


K,y ■■= {x) {x)... ^ (x), <y9 (y), (y),..., ^ (y)) e S{I\ {x, y}). 


We say that a is annular crossing on y? if it satisfies at least one of the 
following three conditions: 

1. There exist distinct elements a,b,c,d G I such that ^ ^ ,d} = («> b, c, d) 

and = (a,c) {b,d). 

2. There exist distinct elements a, b,c,x,y G I such that x and y are in 

distinct orbits of ip and such that ^ = {a,b, c) and ^ ^ i j/} ~ 

(a,b,c) (x,y). 

3. There exist distinct elements a, b, c,d,x,y G I where x and y belong 

to distinct cycles of ip and K,y\^a,h,c 4 } = {a,b,c,d) and = 

(a,c) {b,d) (x,y). 


We say that a is annular noncrossing on ip if it does not satisfy any 
of these conditions. We denote the set of a G 5 (/) which are annular 
noncrossing on ip by 5ann-nc (v?)- 


Lemma 2.13. Let ip^ G S (!) with # (<^+) = I, (p G PM (/) as above, and 
let a G PM (I). Then xipp^a) = 2 if and only if a does not connect the 
eyeles of ip, and a\j is disc noncrossing on ip+. 

Let ip+ G S (!) with ff {ip) = 2, and let a G PM (!) such that ±(/J V a = 
l±i. Then x (Tj «) = 2 i/ and only if a connects one orbit J of ip to exactly 
one other orbit of ip which is not — J (call it K), and a\juj^ is annular 
noncrossing on ip\j\jK- 


See 19 for the proofs. For more on the disc-noncrossing conditions, see, 


e.g. il Bill- The noncrossing conditions for the annulus in the orientable 
case are established in . 

The following lemma can be thought of as the additivity of the Eu¬ 
ler characteristic over disconnected components and the computation of its 
value on connected sums. 


Lemma 2.14. Let vr = {Vi,... ,Vm} G V {L) he a partition of I with m 
blocks, and let a G PM (vr). 

Let ip+ G S {!). If ip+ :< n, then 

x{T,a) = X {t\vi ^ + ■ ■ ■ + X > <^\vj ■ 
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If <y9_|_ has one cycle, then 


= xiv’lvr > + ■■■ + x{v\vrr, ’ ^\vm) “ 2(m- 1) . 


Proof. If y 5 _|_ ^ vr, then every cycle of K (cp, a) is also contained within a 
block of TT. The first part of the result follows. 

Let ifo be a permutation and let t = {a, b) £ S (!) be a transposition 
connoting two cycles of (/pq and two blocks of ±a. Then (a, b) connects two 
cycles of K {ipo, a) (as above), so left multiplying (a, b) and right multiplying 
by (—a, —b) both reduce the number of cycles, that is x {'PoT, a) = x (po, <a)~ 
2. Likewise, if r does not connect cycles of po, but connects blocks of ±a, 
then xiPoT,0() = xiPo^ct) + 2. Thus, if we can construct p+ from pQ by 
right-multiplication by a series of transpositions ti, ... ,Tm £ Sn that do not 
connect blocks of To;, then 

X iP+, ol) = x {po, a) + 2{ff {p+) - # {pq)) . (2) 

Let G S (I) have a single cycle. We note that p'j^ := • • • P+\vm 

differs from only on fc G / where k belongs to a different block of vr 
from (/c), so p'_^pf^ consists of single element cycles as well as cycles 
(ci ,...,Cm) where c* always belongs to a different block of vr than Cj+i. 
We can write such a cycle (ci,..., Cm) = (ci, C 2 ) (c 2 , C 3 ) • • • (cm-i, Cm) as a 
product of transpositions containing elements in different blocks of vr. Then 
since K = p'j^pf^K {p+,a) p-p'p^, the second part of the result 

follows from Q. □ 

2.2 Quaternions 

Definition 2.15. The quaternionsM are an algebra a+bi+cj+dk, a, b,c,d £ 
M, with i,j, k satisfying ij = k, jk = i, and ki = j. 

The real part of a quaternion is 

Re {a + bi + cj + dk) = a. 

The quaternion conjugate is 

a + bi + cj + dk = a — bi — cj — dk. 

The quaternions may be faithfully represented as 2 x 2 complex matrices 

f a + bi c + di 
\ —c + di a — bi 
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We will use 1 and —1 to index the rows and columns. We note that the real 
part is the normalized trace, and the quaternion conjugate is the (complex) 
conjugate transpose of the matrix. 

2.3 Matrices 

We will often use subscripts on matrices, A^. We will often move the sub¬ 
script to a superscript, in brackets: We will use a bracketed superscript 

of 1 to represent the matrix itself, and a bracketed superscript of —1 to rep¬ 
resent the (quaternion) conjugate transpose of the matrix. These notations 
may be combined: 

Definition 2.16. We denote the usual trace by Tr (T) = the 

normalized trace by tr = -Tr. 

Definition 2.17. A bracket diagram is a product in which subexpressions 
may be enclosed in brackets. 

We define the permutation vr G PM {loo) of a bracket diagram on terms 
with subscripts in I by placing a term with subscript oo before and after the 
expression (outside any brackets). For each A: G /qo, let tt {k) be the subscript 
of the term after the one subscripted k, ignoring bracketed intervals. 

Definition 2.18. Let vr, p G PM ([n]^). Then we define a quaternion-valued 
matrix in terms of its entries: 



n 




2:FD(7r)\{oo,dip(oo)}—>-[A^] /cGFD(7r)\{oo} 
/i:FD(7r)\{cX),7r(oo)}^{l, —1} 


We note that there is only one p G PM ([nj^^) corresponding to a given 
FD (p), so p {k) is well defined, even if k is not in FD (p). See for how 
this generalizes the more standard trace along a permutation, as well as the 
proof of the following lemma: 

Lemma 2.19. Consider a bracket diagram on a product of Xk, k ^ I, where 
Re and tr are applied to braeketed intervals. Let it be the permutation of the 
bracket diagram with only the brackets to which Re is applied, and let p be the 
permutation of the bracket diagram with only the brackets to which tr is ap¬ 
plied. Then the value of the bracket diagram is ReFD( 7 r)trFD(p) (-^i) • • • j -^n)- 
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Thus, for example 

Xitr (Re (X 2 *Re (X 3 )) X^Re (X 5 )) Xg = Re^tr^ (Xi, X 2 , X 3 , X 4 , X 5 , Xg) 
where vr = (cx), 1, —4, 6 ) (—2) (3) (5) and (00, 1, 6 ) (—2, 3, —4, 5). 

2.4 Moments and cumulants 

Throughout, (n,X, P) will be a probability triple. 

We will give a definition of cumulants which may be applied to quaternion¬ 
valued random variables. Since the variables do not commute, we wish to 
take the necessary expected values before we perform the multiplications. 

Definition 2.20. The Mobius function (on the poset of partitions) is a 
function /i : V (J)^ —)• C which is given by 




{\{U eTT :U CV}\)1, TT<P 


llyepl 

0 , 


/i(7r,p) : 


otherwise 


The Mobius function has the property that for any n, p € V (/) with 

TT < p: 



See [^, Chapter 10, for more on Mobius functions in general and on 
the computation of this particular Mobius function. 

Definition 2.21. Let Qi,..., Qn : —>• El be random variables. The nth 

mixed moment is the function 

Q-n iQlj ■ ■ ■ ■) Qn} — E iQl ■ ■ ■ Qn) • 

We define a moment associated to each vr = {Vi,..., Vm} G V {n). For 
each 14 £ we define a probability space (11^,74,^^). For each z G F^., 
we define random variables Q( : Hfc —)• BI which are jointly distributed 
identically to the Qj : H —?■ El. The moment associated to vr is the expected 
value on the product of the probability spaces: 



(so the Q'l^ from different blocks are independent). 
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The cumulant associated with n € V (n) is the function 


pS'P(n):p^7r 

The cumulants then satisfy the moment-cumulant formula: 


Q-tt 


fcp. 

pS'P(n):p^7r 


We dehne 


kn ■■= ki^ 


^ ^ /i (vr, 1^) a-j^. 
■KG'P{n) 


If a cumulant has an argument which is independent from all the others, 
then the cumulant will vanish. We note that /c 2 is the covariance of its 
arguments. See for more on quaternionic covariance. 

If a quaternionic random variable is sufficiently symmetric, then its ex¬ 
pected value is real (since quaternionic conjugation is a diffeomorphism 
which preserves volume). All of the quantities we will be dealing with have 
this property: we can see that the expected value of any expression calcu¬ 
lated with Proposition 2.28 will be real. 


Lemma 2.22. Let Qi, ■ ■ ■ ,Qn : II —)• El 6e quaternionic random variables 
with E real for any ii ... ,im £ N• For tt € V {n), let 

C = (oo, 1,... ,n), and let C = Uven Cly Then 

®7r (Ql) • • • ) Qn) — (E O Re)|A/ (Ql, . . . , Qn) 


and, denoting the smallest element ofV by vq: 


k^ {Qi, ...,Qn)= n ^\v\ (^0, C' (^o), • ■ •, (^'o)) • 

vep 


Proof. We use the notation from Definition 2.21 If we have integrated all 


but one of the probability spaces (say, all except for Dfc), then all variables 
(5( with i ^ 14 have been integrated and are therefore real and commute 
with the Q[ with i G 14- Thus if 14 = {4, • • • j im} with 4 < • • • < im, the 
contribution of the Llk is E (Qjj • • • Qim), that is, E (Qi, ■ ■ ■, Qn)) 

(since the expected value is real), giving the first part of the result. 
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Noting that the Mobius function ^ (vr, p) is multiplicative over the blocks 
of p, we may calculate: 

kp {Qli ■ ■ ■ 1 Qn) — 

TT^p V&p 

which gives the second part. □ 


Example 2.23. We calculate that if Qi, Q 2 , Q 3 , Qa satisfy the hypotheses of 
Lemma 12.221 then 


^4 {Qi,Q2, Qs, Q 4 ) 

= E {Q1Q2Q3Q4) - E (Qi) E (Q2Q3Q4) - E (Q2) E (Q1Q3Q4) 

- E (Q3) E {Q1Q2Q4) - E (Q4) E (Qig2Q3) - E {Q1Q2) E {Q3Q4) 

- E {Q1Q3) E (Q2Q4) - E (QiQi) E (Q2Q3) + 2E (Q1Q2) E (Q3) E {Q4) 

+ 2E (Qig3) E (Q2) E (Q4) + 2E (giQ4) E {Q2) E {Q3) 

+ 2E {Q2Q3) E (gg E (g4) + 2e (g2g4) e (gg e {Q3) 

+ 2E (g 3 g 4 ) E (gi) e (g2) - 6 e (gg e (g2) e (g3) e (g 4 ) . 

We note that the real part function is redundant, since the expected values 
of the products are all real. However, it is used in the expression in order 
to give the order of terms in the products. (It would also be possible to use 
the trace function in its place.) 


2.5 Topological Expansion 

The Weingarten function is defined in and the results in the lemma 

that follows are proven there. (See |20[ for more detail on the quaternionic 
case.) The following definition is the quaternionic Weingarten function. 
(Throughout, the Weingarten function will refer to the quaternionic Wein¬ 
garten function. It may be calculated from the real Weingarten function: 
\YgSp(A^) _ -^gO(-2Af)^ ggg fQj. tables of values.) 

Definition 2.24. For tt, p € 1^2 (n), let d (vr, p) := n/2 — ^ (ny p). 

Define matrix with rows and columns indexed by 1^2 {n) by letting Gr (vr, p) 
(_l)”'/2 for 7r,/9 G 'P 2 (n). We dehne the (quaternionic) Wein¬ 

garten function and normalized Weingarten function Wg, wg : T*! —)• M: 

Wg(7r,p) := Gr"^ (7r,p); 

wg (vr, p) := (_ 2 IV)"/ 2 +'^D,p) ^ 
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We note that the Weingarten function depends only on the block struc¬ 
ture of TT V /?. Since all the blocks contain an even number of elements, 
we define an integer partition A (vr V p) whose parts Ai,..., are half the 
sizes of the blocks of vr V p. If A is an integer partition of re, we define 
wg (A) := wg (tt, p) for any vr, p such that A (vr V p) = A. 

For 7ri,7r2 G 1^2 (re) with vri V 7r2 = vr G "P (re) and p,a £ V (re) with 
TT ^ p ^ fj, we define a sort of cumulant of the Weingarten function by 

Cn,p,a-= ^ p(r,f7) Wg(7ri|^, 7r2|y), 

T'.p^T^cr V Gt 

or equivalently, 

wg(7ri|y , 7r2|y) = 

V Gcr T'.p^T^a 

The following is from [^[^, where more details can be found: 

Lemma 2.25 (Collins, Sniady). The Weingarten function may he expressed: 


wg (vr, p) 


^2JYjd(vr,p) ^ 2 JV) 


-(d(cro,o-i)-l- 


k>0 


cro,---,o-keP2(n) 
(JQ —TT —p 


For vr, p G P 2 {n), 


lim wg(vr,p)= IT (-1)* ^ Cj-i + O (A^ ^) 

2 GA( 7 r,p) 


N^oo 


where Ck '■= T^(^k) Catalan number. 

For VTi, vr 2 G V 2 (re) with vri V vr 2 = vr G P (re) and p,cr £ V (re) such that 
TT ^ p ^ fj, the cumulant c,r,p,cr is O (A^ 2 (#(cr)-#(p))^ _ 

Proof. The form for the Weingarten function and its asymptotics may be 
found in [^. 

We define functional S on functions on the positive integers N* —>■ C 
which are eventually zero by S (/) = (—1)^ / {k) and a product g * h 

by 

{g*h){k)= Y. g{\Ui\)h{\U 2 \). 

Ui,U2Qlk]:UiUU2 = lk] 

The product is associative, and S {g * h) = S {g) S (h) (see [^). 
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Let T = G V {n). For m G [t] and for each nonnegative 

integer I, define fm : N —>• C by letting fm {k) be the number paths in 

'P 2 (Ym) (distance metric d) of k steps and total length I (where each step 
is distinct from the previous one). Treating wg as a function of 2N, the 

coefficient of in OveT(^i| v > '^ 2 \y) is 


/i ,••• ,lt :/iH- \-lt=l 

/l,”- - \-lt=l 


which may be interpreted as 


the coefficient of ^ 


in 


{2N) 

k>0 vo,...,Vk€P2(T) 

vo='iri,Vk=Tr2 
^o¥=^l¥=-"¥^^k 


- \-d{v^._i,Vk)) 


(where the Um from the definition of the product * are interpreted as the 
sets on which the km steps of the path on Vm take place). 

The cumulant Cjr,p,(T is the sum over r with p ^ t ^ a oi p {t, a) times the 
above sum. If we consider all occurrences of the term associated with path 
uo,..., Ufc, it will occur in the sum associated with all r with vqV- ■ - Vvk ^ r. 
The sum of the coefficients will be Ylr-p voV---vvk-<T-<cr k- ij-! ^)) which vanishes 
unless pV vqV ■■■ V Vk = (T. Thus the cumulant may be expressed 

’E(-i)‘ E - 

k>0 V0,...,Vk&P2(T) 

V0='!ri,Vl,.=TT2 

pVi;oV---Viifc=(T 


Let uo,..., Ufc be a path in V 2 (n) with uq V • • • V = p of total length 
1. Joining v G V 2 (n) to Vm must reduce the number of blocks by at least as 


much as joining u to p (Lemma 2.3). Thus d (v, Vm) = 'ra/2 — # (u V Vm) > 
# (p) — # (u V p), so inserting it into the path must increase the total length 
of the path by at least 2^ ((p) — # (u V p)). Thus the length of the path from 
TTi to 1^2 must be increased by at least 2 (vr) — (u)) > 2 (# (p) — # (< 7 )). 

The result follows. □ 


The following definition is adapted from in [^|^. See also [^. 
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Definition 2.26. Let Xi,..., Xn ■ ^ —)• Mjstxn (H) be random quaternionic 
matrices. Then the normalized symplectic matrix cumulant is the function 
/ : PM (n) ^ C 

/ (a) = (—wg (ha, (Jvr) 

7rGPM(n) 

X E (ReFD(7r)trFD(7r) (-^1) • • ■ ) ^n)) ■ 
We note that if the matrices have a limit distribution, then 


lim /(a) 

N^OO 

= V' TT Ci-i( lim oEoReotrVwi,...,Xn) 

' ^ \N—^oo j 

7rGPMnc(Q) iGA((5a)V((57r)) 

TT^a 

which is multiplicative in the sense that for r = {Vi ,..., Vm} G V (n), if 
a ^ ±r, then, denoting the matrix cumulant of X,, i G 14 associated with 
a|y, by/fc(a|^J, 

lim f (a) = lim fi(a|,/)--- lim 

The cumulants are multiplicative over independent matrices (see [20] ): 

Lemma 2.27. Let tc : [n] —)> [C] be a word in colours [C], and let Xi,..., X„ : 
n —)• Mj^xN (H) be symplectieally invariant random quaternionic matrices. 
For each c G \C], let the Xk with w{k) = c be independent from the rest of the 
matrices, and have normalized matrix cumculants fc : PM (tc"^ (c)) —>• C. 
Then the normalized matrix cumulants f : PM (n) —)> C vanish unless 
a 4 ±ker (w), in which case 

f («) = /i («L-i(i)) ■■■ fc («L-i(c)) • 

Integrals of traces may be computed in terms of the matrix cumulants. 
See (2^: 

Proposition 2.28. Let Xi,...,X„ : 12 —)• H 6e symplectieally invariant 
matrices with normalized matrix cumulants f : PM ([n]) —)• C. Let Ttr £ 
[n]g^. Let e : [n] —>■ {1, —1}, and let hg G 5 (± [n]) : /c i—)■ e (|A:|) k. Then 


E 


Re tr ( 


— ^_2)X4Re,5£a(5s)-2#(vPRe) jyxiV‘tr,Ssa5s)-2#{iptr) J 

Q;GPM(n) 

where a acts trivially on Too. 
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2.6 Freeness 


Definition 2.29. A noncommutative probability space is a pair (A, </>i) con¬ 
sisting of a unital algebra A and a linear functional 4>i : A —?■ F such that 
(j)i (1^) = 1. (In this paper, F will generally be H, although F = C and 
F = M are more frequently studied.) 

Subalgebras Ai,..., Ac C A are free if (fi (xi • • • Xp) = 0 whenever, for 
all A; G [p], (/>! {xk) = 0 and Xk G A^t^k) where u) : [p] —>• [C] is a word with 
w{l)i^w{2)^ ■■■ ^w(p). 

Sets are free if the algebras they generate are free. (If A is a *-algebra, 
the algebra generated by a set should also be closed under the operation *.) 

Definition 2.30. If we have a set of random N x N matrices 
for arbitrarily large G Z (we will usually suppress the dimension in the 
notation), then we say they have a limit distribution C A, where 

(A, (fi) is a noncommutative probability space, if, for any Ai,..., Ap G A: 

lim E (tr (Aa^ • • • Aa^,)) = (pi (xai • • • xx^) 


and all higher cumulants of normalized traces of products of the Aa vanish 
for large N. 

We may also wish the limit distribution to respect a structure such as 
an involution * : A —>• A, a i—)• o*. We incorporate this into the subscript 
notation: for each A G A, we must have — A G A, and A^ = A_a (where 
here * is the conjugate transpose) and x*^ = x_a. 

In the quaternionic context, we would like the involution to be conjugate- 
linear and reverse the order of multiplication: {qx)* = qx* and {xy)* = y*x*, 
for g G HI and x,y G A. 

Sets of matrices are asymptotically free if their limits in A are free. 


Second-order probability spaces are constructed in [7,14,15 in order to 
model the covariance of traces, in order to study the fluctuations around the 
limits (central limit type results, as opposed to law of large numbers type 
results). 


Definition 2.31. A second-order probability space is a triple (A, (/>i, (^> 2 ) 
where {A, (pi) is a noncommutative probability space and ())2 : A^ —)• F is a 
bilinear function. (In the complex and real cases, (p 2 is required to be tracial 
in each of its arguments.) 

A quaternionic second-order probability space must also be a =i=-algebra, 
where * is an involution A —?■ A, a 1 —)• a*, which is conjugate linear and 
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reverses the order of multiplication. (In the complex case, the involution is 
not required, and in the real case the involution is linear.) 

We say that subalgebras Ai,...,Ac C A are quaternion second-order 
free if they are free; and letting n : [p] —)• [C] be a word with either u (1) 7 ^ 
u ( 2 ) / • • • / u (p) 7 ^ u ( 1 ) or p = 1 and u : [g] —)■ [C] be a word with either 
u ( 1 ) 7 ^ u (2) 7 ^ • • • 7 ^ u (g) 7 ^ u ( 1 ) or (7 = 1 , when p 7 ^gorifp = (7 = l and 
u(l) 7 ^ u(l): 

4>2 (a;i • • •Xp,pi • • - pq) = 0 

and when p = q > 2: 


4>2 (xi • • • Xp, Pl • • • Pq) 

P ^ P-1 P 

= Re {cl)i (xiPp+i-i)) + ^ X] n (*^1 (xiPp+i-k-i)) 

i=l k=li=l 


p . P-1 P 

- 2 (xiv*)) + {xiVk+i)) ■ 

i=l k=li=l 


(All indices are taken modulo the appropriate value to fall in the appropriate 
interval.) 

We note that if the subalgebras are quaternion second-order free, then if 
Xi and Uj do not come from the same subalgebra, (fi {xipj) = (fi {xiy*^ = 0 
(since they are free) so any term containing such an expression will vanish. 

Sets are quaternion second-order free if the ^-algebras they generate are 
quaternion second-order free. 


The expression for the value of 4>2 is motivated by Proposition 4.4 For 
the expression in the complex case, see [7, 14 15 , and for the expression in 
the real case, see [^. Like the complex and real cases, the terms correspond 
to spoke diagrams (see Figure [TG] ), but in the quaternionic case they appear 
with different weights and signs. 


Definition 2.32. We say that matrices have a second-order limit 

distribution {xA};ygA ^ ^ if this set is a limit distribution for {A,(fi), and 
in addition 


(Tr (X,. ... , Tr ■ ■ ■ A\J) 

</’2 (xAi ■ ■ ■ ^Ap) S^Ap^x ■ ■ ■ 2^Ap_|.q) ) 

and all higher cumulants of unnormalized traces vanish for large N. 
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Figure 2; The four types of spoke diagrams which appear in the quaternionic 
case. The solid polygons represent the cycles of tr, and the dotted lines the 
cycle of Re. The upper left is a sphere {k = p where k is the signed offset 

X = 2), the upper right a torus (A: = 2, y = 0), the 


defined in Lemma 4.3 


lower left a projective plane {k = —1, x = 1); and the lower right a Klein 
bottle {k = —2, X = 0)- 
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In the quaternionic case, we require that A has a conjugate-linear invo¬ 
lution *: A —)■ A, X I—)■ X*, and for each A G A, —A G A, and X_x = and 
X-X = x\. 

We say that matrices are asymptotically quaternion second-order free if 
they have a second-order limit distribution and their images in the limit 
distribution are second-order free. 

We note that if matrices have a limit distribution, then if we expand a 
product of traces in cumulants: 

Jim Ott (tr (Xi),..., tr (X^)) = Jim XU k\v\ (tr(Xq),...,tr . 

p^TT Vdp 

Since all cumulants km with m > 1 vanish as X —)■ oo, the only remaining 
term is 

lim E(tr(Xi))---E(tr(X„)). 

N^oo 

If they have a second-order limit distribution, then the difference from this 
value is O (X“^). 


3 Cumulants 

In this section we show a sense in which the cumulants of traces of products 
of symplectically invariant random matrices are sums over connected terms, 
and look at bounds on the order of cumulants of traces of expressions in 
independent matrices with limit distributions. 

Proposition 3.1. Let 0 < ni < • • • < Um, and let Xi,...,X„^ : II —)■ 
MmxN (H) symplectically invariant random quaternionic matrices. Let 
w : [rim] —)■ [C] he a word in colours [C], and let the matrices Xk with 
w (k) = c be independent from the other matrices, with matrix cumulants 
fc-.PM {w-^{c))^C. 

Let 

c = (oo, 1 ,. . .,nm) 

and let 

ip = (oo) (1 ,..., ni) • • • {Um-l + 1, • • •, rim) ■ 

Then, denoting the cycles o/FD (tt) by tti, , Tip, and for a {p), letting 


21 



^ = {Ufcel/^fc : V e ^}- 


kn (tr (Xi • • • ^m),. ■ ■, tr • • • ^um)) 

— (•_ 2 )X(C>“)- 2 +x(“>'^)-#(«) jyxiv,o:)-2m+x{a,-n-)-#{a) 

oSPM(n),7rSPM(n) 
p:(<5a)V(57r):;<p:;<it ker(ti>) 
o'SP(p):o':^it ker(ui) 

±pVQ:VpV<3-=l±[„^] 

^C(^Q,)V(5jr),((5Q:)V((57r),p^(T O tr)^^ (^1) • • • ) ^rim) j • • • ) ° ^'^)'7rp (^1) • • • ) ■ 

Proof. For t eV (m), we define r e P (rim) by f := {UfceP + 1; ^fc] ^ ^ ^ 'ri¬ 
nsing Proposition 2.28[ we compute the summand in the cumulant corre¬ 
sponding to r: 

^(r,l^) (_2)X(C,«)-2^x(v^,a)-2m 

QGPM(fAker(ti;)) 


7rGPM(fAker(ip)) 


n wg((5Q;|^f;, 

k [/ETAker(ii;) 


\ PGTAker(tt;) 


(The Euler characteristics are calculated using Lemma 2.14 Ylv&r X {Pi Q^) = 

X {if, a) and J^ver (x (Cly , «) - 2) = X (C, «) - 2.) 

We collect terms associated to a fixed a G PM (ker (re)) (there will be 
one for each t € V (m) such that ±f ^ a). We expand the Weingarten 
function and expected value in cumulants: 

Wg ((5 q!|^jj , (JvrAjy) ^ ^ ^(Sa)V(STr),(Sa)V(S7T),p 

U G'rAker(ti;) p:((5a)V(57r):Kp^difAker(ip) 

(so all p ^ ± ker (tc), and a term associated with p appears in the expansion 
for r if ±f ^ p). We likewise expand the expected value in cumulants: 


n IE ° t^)FD(7r|j^^) (^ 1 ’• • • ’ 

yGrAker(it;) 


^ k. ((Re o tr)^^ (Xi,..., ,..., (Re o tr)^^ (Xi, ...,Xnj) 


cr&'P{p) 
itCT:;< itf Aker(-u)) 
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(so all a have a ^ ±ker (w) and the term associated with a appears in the 
expansion for r if ±r ^ a). Collecting terms associated to a hxed a, p, and 
a, the sum of the Mobius functions (noting that p (r, 1^) = p l±[nm])) 
is 

r:±'f^±9?VQ:VpV(T 

which vanishes unless ±(p V a V pV a = □ 


Corollary 3.2. If independent sets of random matrices have limit distri¬ 
butions, the the *-algebra they generate also has a limit distribution, and if 
they have second-order limit distributions, then the *-algebra they generate 
also has a second-order limit distribution. 

If m > 1, then the mth eumulant is 0{N~'^), and in the notation of 
Proposition 3.1. terms of order have a G PM^c (<p) and n G PMnc (a)- 


Proof. By Lemma 2.9 we know that y {(p, a) < 2ff (±(p V a) and y {a, tt) < 


2ff (±a V vr), so y (tp, Q:)+y (a, tt) < 2ff (±(p V a V t:)-\-2^ (i«) (Lemma |2.3[ ). 
Similarly, {5a) V ((57r) ^ zbtp V a V tt, so y ((p, a) + y (a, vr) + 2^ (p) — 
2ff (((^a) V (<5vr)) < 2ff (±(p V a V vr V p) + 2^ (±a). Thus if m = 1, the 
limit exists, so the matrices have a limit distribution 

Since the X^, w {k) = c have a second-order limit distribution, then if 
block P G cr has |P| > 1, the eumulant associated with the block V is 
(9(jV“l'^l), and limAr^.oo vanishes for |P| > 2. (Mixed cumulants 

always vanish, since the terms are independent.) In a nonzero term, ±(p V 
a V p V d = so if ±(p V a V vr V p has more than one block, each 

block must contain at least one vr^ contained in a block of & containing 
more than one element. We calculate that if m > 1, the mth eumulant is 
O {N~'^). Furthermore, if m > 2, then either a block of a with more than 
one element has more than two elements, in which case the limit of iV™ times 
the eumulant vanishes, or the blocks of cr with more than one element all 
have two elements, in which case more than ((p V a V vr V p) elements are 
contained in such blocks, and the eumulant is O ^ for some m' > m. 

Thus the generated algebra has a second-order limit distribution. 

A term from the eumulant which does not vanish faster than (or 

when m = 1) must have y (tp, a) = 2ff (<p V a) and y (a, vr) = 2# (a, vr), 
which demonstrates the second part. □ 


Example 3.3. The diagram shown in Figure[^corresponds to an a G PMnc (<p)j 
as would appear in the expansion of 


fes (tr (X 1 A 2 X 3 X 4 X 5 X 6 ), tr (A 7 X 8 A 9 A 10 X 11 A 12 X 13 ), tr (XuXisXm)). 
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Since ^ip V a = {± [13] ,± [14,16]}, either vr ^ a or tt connects at least 
one cycle of a from ± [13] to one cycle from from ± [14,16], as well as their 
paired cycles (so y (a, vr) = # (a) — 2). In the first case, either p or a must 
connect ± [13] and ± [14,16]. 

We calculate 

K{C+,a) 

= (1,16,15,13, -3, 6 ) (- 6 , 3, -13, -15, -16, -1) (2,4, -10) (10, -4, -2) 

(5) (-5) (7) (-7) (8) (-8) (9) (-9) (11,12) (-12, -11) (14) (-14), 

so y (C, a) = 0 , contributing a factor of 

4 Second-Order Freeness 

Lemma 4.1. Let p > 1, let ifj^ = (1,... ,p) G Sp, let tp- = (—1,..., —p), 
and let p = p+pZ^ ■ Let a G PMnc (p) have no cycle containing both k and 
k + 1, k G \p — 1\. Then a has at least one cycle with only one element. If, 
in addition, p is not in the same cycle of a as 1, then a must have at least 
two cycles with only one element. 

Proof. We can find an a G \p] such that a~^ (a) > a (if all elements 
of Ip] are all fixed points, the result follows, and if a~^ (a) < a, then 
a~^ (a~^ (a)) > a~^ (a) since a is disc standard). Then a~^ (a) / p (a), so 
there are p (a ),... ,p^ (a) which do not share an orbit with a or any other 
elements of \p] (by the disc-nonstandard and disc-noncrossing conditions, 
respectively). If any of these is not a fixed point, then we may repeat the 
argument with a, k' < k. Thus there must be at least one fixed point. 

If p and 1 are not in the same orbit of a, then we may repeat the 
argument taking p (a) as a, which gives us a distinct set of elements, so we 
may find at least one more fixed point. □ 

Proposition 4.2. Independent sets of symplectically invariant quaternionic 
random matrices with a limit distribution are asymptotically free. 

Proof. Let Ai,..., Ac be independent sets of symplectically invariant quater¬ 
nionic random matrices, and let Xi,..., Xp ; 11 —)■ El have E (tr (W^)) = 0 
and Xk G where re : [p] —)• [C] is a word where tc (1) / rc (2) / • • • / 

w (p). We consider the value of 

lim tr (Xi ■ ■ ■ Xp). 

N^OO 
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According to Proposition [3^ (using the same notation) when m = 1, surviv¬ 
ing terms will have a G PM {ip) and -k < a (since x («; tt) < 2# (±a V tt)). 


By Lemma 4.1, there is a cycle of a, and hence of tt, which contains only 
one element, so a term E (Re (tr (X^))) = 0 appears for some A; G [p]. Thus 
the expression vanishes. □ 


Lemma 4.3. Let 


p+ = (l,...,p)(p+l,...,P + g) G Sp+q. 


Let p- = (-1,...,-p) (-(p -t- 1)(p-I-g)) and p = p+P- ■ Let 
a G PMjic {p+P~l^') have no fixed points and no cycles containing neigh¬ 
bours a, p+ (a). Then p = q, and a is a pairing, with pairs of the form 
(i, (p“* (± (2p -I- 1) — k)), where k G ± [p -t- 1, 2p] can he thought of as a 
signed offset (see Figure^2^. 

// = (1,..., 2p) G S 2 p, then x (C+)«) depends on k as follows: 


x(C+,a) = 


2 , k = p 
1, k = -l 
0 , otherwise 


Proof. Since a G PMnc {p+pZ.^), then for some choice of sign, it does not 
connect I := [p] U (± [p -|- l,p -|- g]) to —I. Throughout the rest of this proof, 
we will implicitly consider p and a on 1. Then a is annular noncrossing on 


p as described in Definition 2.12 


We show first that a cannot contain any a and b in the same orbit of 
p. We know there is at least one cycle of a which contains an x, y G I in 
distinct cycles of p, and G S'disc-nc (-^a;,!/) (the second and third an¬ 

nular crossing conditions imply the first and second disc crossing conditions, 
respectively, for cycles not containing x and y, and the annular nonstandard 
conditions and the first annular crossing condition imply the first and second 
disc crossing conditions, respectively, when one of the cycles involved is the 


one containing x and y). Then, following the proof of Lemma 4.1 if a cycle 


of a contains two elements from the same cycle of p, then a cycle of ot must 
contain either a singlet or two neighbours in \x,y 

Thus a must be a pairing on I, in which each pair contains an element 
of [p] and an element of ± [p -|- l,p -|- g] (so p = q). Define a signed offset k G 
± [p] by letting ± (2p -|- 1) — A: be the partner of p. Then for f G \p], i is part¬ 
nered with p~^ (± (2p -|- 1) — A:). (Otherwise, let i be the first element of [p] 
not partnered with (A;); the partner of i will be (± (2p -|- 1) — A:) 
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for some 0 < m < p. Letting x = i — 1 and y = (± {2p + 1) — k), 

the third annular noncrossing condition with Xx^y implies that the elements 
in {</:>“* (± {2p + 1) — /c),..., (i {2p + 1) — /c)} may not share cy¬ 

cles with elements outside that set. Since they are all elements of the same 
cycle of (f, they may also not share cycles with each other, a contradiction.) 

We may calculate K (C+,«) for several cases. If k = p, then 

FD {K (C+, a)) = (1, 2p - 1) (2, 2p - 2) • • • (p - l,p + 1) (p) {2p) 
so x (C+, a) = 2. If /c G [p — 1], then 
FD (iF(C+,a)) 

= (1, 2p — A: — 1) (2, 2p — A: — 2) • • • (p — /c — l,p -|- 1) (p — /c,p, 2p — k, 2p) 

{p — k + l,2p — 1) ■ ■ ■ {p — l,2p — k + f) 


so X (C+) a) = 0. If A: = —1, then 

FD {K (C+, a)) = (1, - (p + 1)) (2, - (p -h 2)) • • • (p - 1, - (2p - 1)) (p, -2p) 
so X (C+) Oi) = 1. If A: G [—p, —2], then 
FD {K{C+,a)) 

= (1, -{2p + k + 2)) (2, - (2p + A: + 3)) • • • {-k, - (2p - 1)) 

{-k + 1, -2p, 2p + k + l, -p) {-k + 2, - {p + 1)) ■ ■ ■ {p - 1, - (2p -h k)) 

sox(C+,a) = 0. □ 

Proposition 4.4. Let Ai,..., Ac be independent sets of symplectically in¬ 
variant quaternionic matrices with second-order limit distributions. 

Let Xi,...Xp,Yi^... ,Yq : Ll —)■ Matxtv (H) be random quaternionic 
matrices with E(tr(Xfc)) = E(tr(yfc)) = 0 for all k, and such that Xk G 
^u{k) and Ffc G for words rt : [p] —)■ [C] with either u(l) / «(2) / 

• • • / rt (p) / tt (1) or p = 1 and r : [g] —)• [C\ with either v {!) ^ v (2) / 

• • • / r (g) 7 ^ r (1) or q = 1. 

Then if p / q, or if p = q = 1 and rt (1) / r (1); 

lim k 2 (tr (Xi • • • Xp ), tr (W • • • P,)) = 0 


26 



and if p = q > 1: 


lim k 2 {ti{Xi---Xp),tT{Yi---Yg)) 

N^oo 

P . P-1 P 

= TT lim E (Re (tr {XiYp+i_i))) + 7 TT lim E (Re (tr {XiYp+i_k-i))) 

i=l k=l i=l 

- 2 n +4 E n 

i=l k=li=l 

where indices are taken modulo p. 

Remark 4.5. We note that if u (k) 7 ^ v (1), then 


E (Re (tr = /„(fc) ((fe) {-k)) {{1) (-/)) 

= E (Re (tr (X^))) E (Re (tr (Yi))) = 0, 

so any summand containing a term which matches matrices with u (k) 7 ^ v (1) 
vanishes. 


Proof. Let 


let 


ptr ■= (oo) ,p) {p+l,...,p + q) 


(fRe ■■= {oo,l,...,p + q), 
and let u) : [p + (?] —)■ [C] be given by 


w {k) = 


u (k) 


k £ Ip] 


V {k — p) k£\p + l,p + q] 


Up = q = 1 and tt (1) 7 I: u (1), then k 2 (Tr (Xi) ,Tr (Ti)) is the covariance 
of independent quantities and hence vanishes. 

U p ^ q or p = q > 2, then applying Propositionto the left-hand side 
of the desired covariance in the statement of the proposition, we note that if 
a does not connect the blocks of ±ip, then at least one of ajj^ and a|[p+i p+gj 
must contain at least two single-element cycles (Lemma 4.1). If tt connects 


the blocks of ±(/?, then in an asymptotically nonvanishing term tt must con¬ 
nect at most two blocks of ± 0 ; (since y (a, vr) < 2ff (±a V vr)), and u = Op. 
Thus a term of the form E (Re (tr (X^))) = 0 or E (Re (tr {Yk-p))) = 0 must 
appear in the expansion, where k is the element from a single-element cycle 
of a which vr does not connect to another cycle, and so the term vanishes. 
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Likewise, if /? or cr connects the blocks of ±(/?, then vr must not connect any 
cycles of a, and there may not be more than one block of a containing more 
than one element, so there must again be a term E (Re (tr (X^))) = 0 or 
E (Re (tr (Yfc_p))) = 0 where k is the element from a single-element cycle of 
TT which appears in a single-element block of a. Thus all terms in which a 
does not connect the blocks of vanish asymptotically. 

In an asymptotically nonvanishing term in which a connects the blocks 
of ±(/9, vr does not connect cycles of a, p = (5a) V (dn), and a = Op. We 
calculate: 

hm k 2 {Tv{Xi---Xp),TT{Yi---Yg)) 

N^OO 

oSPMnclv’tr) 7re5nc(a):7r^Q; yieA(((5o)V(<57r)) 

±<PVo=l±[p+5] 

X (E o Re o tr)pj-)j-^^ (^ij • • • j ^pi hi, • • •, Yg). 

If a contains any single-element cycle, then E (Re (tr (X^))) = 0 or 
E (Re (tr (Yfc))) = 0 must appear for some k. Thus we may restrict our 
attention to a which satisfy the hypotheses of Lemma |4.3[ 

Since such a are pairings, any tt < a has cycles of length 1 or 2. If 
any cycle of vr has length 1 the term again vanishes. Thus vr = a for any 
nonvanishing term. 

For such a term, 5a = 5 tt, so A {{5a) V (^tt)) = (1,..., 1). Thus the 
contribution of the product of Catalan numbers is 1. Lemma |4.3| gives the 
possible values of x((/9tr,a), and hence the power of —2, for the possible 
cases of a. □ 


This culminates in the following: 


Theorem 4.6. Independent sets of symplectically invariant random quater- 
nionic matrices are asymptotically quaternionic second-order free. 


Proof. Asymptotic freeness is shown in Proposition 4.2 second-order limit 


distribution follows from Corollary 3.2, and the formula for the asymptotic 


covariance of traces is given in Proposition 4.4 


□ 


5 Specific matrix models 

We demonstrate in this section that four important matrix models have 
second-order limit distributions, so the results of this paper apply. We 
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consider upper bounds on cumulants of traces, which are of the form: 


km. itr [X 


d£(i)) 




■ • '■») . ■ ■ ■. ) 

aePM(n,rt) 

X Yl f {a\v) ■ (3) 


7r£'P{m):7ry(y 


Veii 


We also give some formulas for the contribution of the four matrix mod¬ 
els to the asymptotic covariance. We may construct a centred element of 
an algebra of random matrices from any element by subtracting the ex¬ 
pected value of its trace: X — E (tr (X)). Then the factors appearing in the 
summands for the asymptotic covariance are of the form 


(E o Re o tr) ([X - E (tr (X))] [X - E (tr (X))]) 

= E (Re (tr (XX))) - E (Re (tr (X))) E (Re (tr (X))). 

(The value of the function E o Re o tr is unchanged when the terms in the 
product are cycled, and E commutes with Re.) Since matrix cumulants are 
asymptotically multiplictive, this is the sum over cumulants associated to 
dia gram s which connect the intervals corresponding to X and X (see Fig- 
- x(^(i)) ... and X(±i) = where 


ure 


5.1). IfX = Z) 


_ -‘m+l ^m+n 

Zi,. ■ ■, Zm+n have matrix cumulants / : PM (m -|- n) —)> C, then denoting 
C = (1,..., m -|- n) (— (m -|- n),..., — 1 ) and vr = {[m], [m -|- 1 , m -|- n]}: 


lim (E (Re (tr (XX))) - E (Re (tr (X))) E (Re (tr (X)))) 


E 

oePMnc(C):±’rV«=l±[„,+„] 


lim / {6ea6e) ■ 

N—^OO 


Example 5.1. If, in the calculation of lim^r_,.oo k 2 (Xi • • • Xp, Xi • • • Yp), we 
encounter a term X^ = IXi • • • We, paired with a term 1 )* where T) = 
1 X 7 • • • VFio, we consider connected, noncrossing diagrams such as the one 
in Figure 5.1 which we can think of as an enlargement of one of the spokes 


in the diagrams in Figure 2.6 


The diagram in Figure [TT] contributes 


Re (tr {DiD^DeDlo)) Re (tr {D 2 )) Re (tr {0^0^)) Re (tr {D*^D*gD*s )). 
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Figure 3: An enlargement of a term from a spoke between terms = 
fFi • • • VFe and Y; = • • • VFf. 


Throughout, a standard quaternionic Gaussian random variable is a ran¬ 
dom variable + ^ 2 j + where the are N (O, |) random variables 

and G : n —)■ MmxN (H) is a random M x N matrix whose entries are i.i.d. 
standard quaternionic Gaussian random variables. The matrix cumulants 
are calculated in 


20 


5.1 Quaternionic Ginibre matrices 

Ginibre matrices are hrst defined in [10[ . 

Definition 5.2. A quaternionic Ginibre matrix is a random matrix Z = 
-^G with M = N. 

VN 

The matrix cumulants of Ginibre matrices are 


J 1, a G 1^2 (± N), sgn {a (k)) = —sgn {k) for all k 
\ 0, otherwise 


Since the cumulants are multiplicative, 


f {a\v) 

7v^V{m):'kya V^Gtt 


1, ±ipV a = l±n^ 
0, otherwise 


so in any non-vanishing term of ([^, x (y?, a) < 2. Thus the *-algebra gen¬ 
erated by a Ginibre matrix has a second-order limit distribution. 

A noncrossing pairing must connect elements with oppositive values of 
e. The number of such noncrossing pairings depends on the values of e. 


5.2 Gaussian symplectic ensemble matrices 

Definition 5.3. A Gaussian symplectic ensemble matrix, or GSE matrix, 
is a matrix {G + G*) with M = N. 
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The matrix cumulants of GSE matrices are 



1, a £V 2 (n) 
0, otherwise 


Like Ginibre matrices, cumulants are multiplicative, so GSE matrices 
have a limit distribution. 

Noncrossing pairings on the disc are enumerated by the Catalan num¬ 
bers. If X = T™ - E (Re (tr (T™))) and E = - E (Re (tr (T^))), then the 

contribution of the spoke is 



m, n even 
m, n odd 
otherwise 


lim E(Re(tr(XE))) 


N^OO 


5.3 Quaternionic Wishart matrices 

Wishart matrices are first defined in 25 . 

Definition 5.4. A quaternionic Wishart matrix is a matrix W = ^G*DG, 
where D G MmxM (H) is a deterministic matrix. We may also consider 
a set of Wishart matrices with different deterministic matrices but the 
same (not independent) G\ Wi = ^G*DiG ,..., Wn = ^G*DnG. 

The matrix cumulant of Wi ,..., Wn is 


/ (a) = (Re o tr)pj 3 (^_i) (Di,..., . 


As with Ginibre and GSE matrices, the cumulant is multiplicative, so 
Wishart matrices have a second-order limit distribution. 

The value of an expression with Wishart matrices depends on the val¬ 
ues of the Dk- If Dk = Im for all k, then it may be expressed in terms 
of the number of noncrossing hypermaps, which are counted by Catalan 
numbers. In this case, if X = Wi ■ ■ ■ Wm — E (tr (Wi ■ ■ ■ Wm)) and Y = 
Wm+i • • • khm+n IL (tr {Wm-\-i''' Wm-i-n))) then 


E (Re (tr (XT))) = - GmCn- 


5.4 Haar-distributed symplectic matrices 

Definition 5.5. A symplectic matrix is a matrix U G Mj^xN (H) such that 
UU* = In- The symplectic matrices Sp (X) form a compact group, which 
has a Haar measure of finite measure. When this measure is normalized 
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to total measure 1, it is a probability measure. A Haar-distributed ran¬ 
dom symplectic matrix is a matrix distributed according to this probability 
measure. 

The matrix cumulants of Haar-distributed symplectic matrices are 

f(a) = l (")))’ sgn(a(A:)) =-sgn(A:) 

^ 1 0; otherwise 

As calculated in Lemma 12.251 

/i (vr, Im) Il f {a\y) = 0 
7rGP(m):7r^CK V 

which, combined with the exponent on N in (© (where < 

2ff V a)) gives us that the mth cumulant of unnormalized traces is 
SO Haar-distributed symplectic matrices have a second-order 
limit distribution. 

Since U* = 17“^, any product of U and U* is of the form U"', n G Z. If 
X = — E (tr (H™)) and Y = U"^ — E (tr (H"')), then since a noncrossing 

a must take any k to an a{k) such that £{a{k)) = —£{k) (i.e. from the 
other term), m = —n and: 

lim E(Re(tr(Ay))) = 

7V->-oo 
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